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Abstract 

The Backlund transformation (BT) for the Camassa-Holm (CH) equation is 
presented and discussed. Unlike the vast majority of BTs studied in the past, for 
CH the transformation acts on both the dependent and (one of) the independent 
variables. Superposition principles are given for the action of double BTs on the 
variables of the CH and the potential CH equations. Applications of the BT and 
its superposition principles are presented, specifically the construction of travelling 
wave solutions, a new method to construct multi-soliton, multi-cuspon and soliton- 
cuspon solutions, and a derivation of generating functions for the local symmetries 
and conservation laws of the CH hierarchy. 


1 Introduction 

The original Backlund transformation (BT) arose in the context of differential geometry 
of surfaces in the 1880s [3]. In the modern era, BTs have been recognized as playing 
a central role in the theory of integrable differential equations [3H |59l |58]. Their pri¬ 
mary application is as a method to generate explicit solutions, exploiting the so-called 
superposition principle, an algebraic rule to “combine” two solutions obtained by BTs 
(from a given initial solution). However, in recent work [52] we have also shown how to 
derive local symmetries and conservation laws directly from BTs. There is also a deep 
relationship between BTs and the associated linear systems of integrable equations. 

The Camassa-Holm (CH) equation [9l[T0] is by now recognized as one of the archetypes 
of integrable equations. It has (weak) “peakon” solutions — solitary waves with discon¬ 
tinuous hrst derivative at their crest — and numerous other types of travelling wave solu¬ 
tion, including solitons (smooth solitary waves), cuspons and various periodic structures 
gni miiziiEi 1331 sa EEi Ea EH El . The integrability of the CH equation was already hrmly 
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established in [9], where a Lax pair and a bihamiltonian structure were given, and much 
further evidence for this has accumulated since then. There is an inverse scattering for¬ 
malism [nun], explicit formulas can be found for multipeakon,multisohton, multicuspon 

and soliton-cuspon solutions [BDl IH El El 1221 ED ED SHI SD ESI SSI ESI SD ED E21SD ES] 
there are an inhnite number of local conservation laws [2S1 ED ED [2S1123 [2D ESI ED EDI 12S] , 
and there is a rich algebra of symmetries [SSIED12S1125112D ES]. Other signihcant works 
on CH include studies of the stability of peakon and other exact solutions ED ED ED ED 
and interesting numerical studies [321 SSI ED E2]. 

The aim of this paper is to fully explore the theory of the BT for the CH equation. 
In [HD], one of us constructed a BT for the associated CH (aCH) equation, an equation 
related by a (held dependent) change of coordinates to the CH equation, and used this to 
construct some solutions of CH which could be regarded as superpositions of 2 travelling 
waves. However, this work was incomplete; an integration was required to reconstruct a 
solution of CH from a solution of aCH, which, in general, could not be done explicitly, 
severly limiting applicability. In the current paper we resolve this and other problems. 
The BT of CH differs from standard ones (for example, those of KdV and Sine-Gordon) 
in that it involves a transformation of both the dependent and one of the independent 
variables. However, remarkably, there is a nonlinear superposition principle for both of 
these transformations, which we develop and apply to the generation of multisoliton, 
multicuspon and soliton-cuspon solutions, as well as to the derivation of symmetries and 
conservation laws for CH. The action of the BT on both dependent and independent 
variables is not unique to CH; a similar situation exists for the Dym equation, which also 
exhibits nonanalytic solitons [HDIHD. 

The structure of this paper is as follows: In section 2 we recap the known results for the 
aCH equation. In section 3 we use them to derive the BT for CH. Section 4 discusses the 
various forms of superposition principle. In section 5 we use the BT to obtain travelling 
wave solutions. The BT is used to construct soliton and cuspon solutions from which 
the standard peakon solutions can be obtained in a certain limit. Alas it does not seem 
to give a direct construction of peakons. However, various other unphysical solutions 
are also obtained. In section 6 we use the superposition principle to obtain cuspon- 
cuspon, soliton-soliton and cuspon-soliton solutions. In section 7, following [53], we use 
the BT to construct the conservation laws and symmetries of CH. Section 8 contains 
some concluding remarks. 

2 Previous results 

The Camassa-Holm equation (CH) [9] is 

rut + 2uxm + urux = 0, m = u — Uxx, (1) 

or equivalently 

Ut — Utxx + “inUx — UUxxx — 2.UxUxx = 0. (2) 

By translating u and performing a Galilean transformation x —?> a; — ct it is possible to 
introduce linear transport and linear dispersion terms into the equation, see for example 
|2Uj . All the results we present here can be generalized for the full class of equations 
considered in [20] . 
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Writing u = Vx and integrating once, we obtain the potential Camassa-Holm eqnation 
(pCH) 

So 1 2 ^ 

'^txx "h 2^* '^x'^xxx 2^3; 0 ! 

or, eqnivalently, 

+ ]^vl + VxUx - - n)^ = 0 , n = v - Vxx ■ 


Evidently n is a potential for m, m = Hx- 

In [60] eqnation ([T]), nnder the assnmption m > 0, was transformed to the associated 
Camassa-Holm eqnation (aCH) 


2pr = —p^U^ , U = 


P Pr 


2 \ p 


+ P , 


with the help of transformation 


p = y/m , = -pdx — -pudt , r = t. 


This transformation implies 

dx 2 dx 

p ’ dr 

A BT for aCH was fonnd in 1601: 


= u 


m 


= 0 , 




= 1. 


p ^ p — 


u ^ u 


2sx 

p{p - s^) ’ 


where s satishes 


5? 


St 


a 

I h p , 

pa p 

2 Px / \ 

— s -I- s -|- a(a + u) . 

P 


The following nonlinear snperposition principle was also given: 


p ^ p — 


/ (q - f3){af3 - SgSfs) 

V /3s„ - 




(3) 

(4) 

(5) 

( 6 ) 

(7) 

( 8 ) 


where Sj^ are the solntions of fl6H7|l with parameters a and (3 respectively. 

In [55] the BT was used to hnd an inhnite number of symmetries for aCH. These are 
given by the generating symmetry X = where 


p{s^a^ + 


Q^ = 


(2)^ 


Q:(Sa^ - Sa’) 


(2)^ 


Q“ = 


2Sa ^ -h 2Sa^ + PU^ 

Ab _ A2) • 

'-’Q *q: 


(9) 


Here Sa \ Sa ^ are two different solutions of (16H7|) for the same parameter a. This symmetry 
depends upon a; expansion in a (formal) power series in a gives the inhnite hierarchy of 
symmetries. 
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3 The Backlund transformation for the Camassa-Holm 
equation 


In this section we obtain the BT for CH and pCH from the BT for aCH. With the help 
of (jl]) we write the BT ([5]) , ([6]) , ([7]) as 


u ^ u — 2a — 


2a{uxS — 

— a 


au) 

2 


( 10 ) 


where s satisfies 


Sx — 




-{m + a) , 


u 

2a 


— UxS + 2 


um) . 


(11) 

( 12 ) 


This system for s is eqnivalent to the Lax pair for CH. Note ffT2|) can be simplified with 
the help of (ITT]) and ([I]) to 


St — T 2aSx — u^x — • 


(13) 


In light of (jl|) the BT for CH must also involve the independent variable x. Using the 
first equation in (Il|), the change of the independent variable is 


Xy 


X = 


Pnew P, 

- - Me 

2s^ 


P{P - s^) 
2ds 




= In 


a 


a 


s + a 


+ /M 


In moving from the third to the fourth line here the formula for in (jH]) is used in the 
denominator but not in the numerator. The integration leaves undetermined an arbitrary 
function /(r). Using the second equation in (j4]) it is straightforward to show this must 
be a constant, which can be taken, without loss of generality, to be zero. Thus the effect 
of the BT on the independent coordinates is 


X —)■ a: + In 


s — a 
s + a 


t —y t. 


(14) 


There is no guarantee that this mapping will be a bijection. We will see later an example 
in which the BT generates several solutions out of one, in the case that this mapping is 
not 1 to 1. 

Using (j5|) and (j6|) it is straightforward to write down the BT for the field p 

p^- -^ (15) 

ap 
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and hence also for the held m = = u — Uxx 


m —)■ 


(s^ — 
a^m 


2 


(16) 


Further calculations give the action of the BT for the pCH helds v (satisfying u = Vx) 
and n = V — Vxx'- 


V ^ V + 


n —)■ n — 2 s , 
2a{aUx — us) 




(17) 

(18) 


As mentioned above, the BT can be generalized for the full family of equations from 


C-lUx T C2Uxxx T C^(Ut T ^UUx) C^iUtxx T UUxxx T ‘2UxUxx) ) 

where ci, 02 , 03,04 are constants. (This generalized equation is referred to in [26l El] as 
the “CH-r equation”.) The BT is 


SyTq — CKyTs 
syTJ + ofyTa 

Here s satishes 



u ^ u — 2a — 


2c^a^u — 2c^asUx + oick^ + 02S^ 
03 — 045 ^ 


^x 


St 



a 2c^u — 2 c 4 Uxx + 2 C 3 CK + Oi 
2 2040 — 02 


^Uxx T 2 q^S,j; — USx — UxS. 


( 20 ) 

( 21 ) 


Equation flT^ includes the KdV, CH, and Hunter-Saxton (HS) [29] equations. The KdV 
equation can be obtained by putting oi = 04 = 0. The HS equation 

Utx T 2 ^^ ^^xx 0 (22) 


can be obtained by putting Oi = 02 = 03 = 0 and integrating with respect to x. The BT 
in this case is 


2 ck 2aUx 

X ^ X -, u ^ u -2q!, 


where 


Sx 

St 


Ur, 


2a 2 ’ 

dUxx + 2q;S3; — USx — UxS. 


(23) 

(24) 


4 The double Backlund transformation and super¬ 
position principles 

In this section we discuss double BTs for CH and pCH. We also show the superposition 
principles for these equations. 

As we saw in the previous section, a BT (which acts on the CH helds M,m,p = ^/m, 
the pCH helds v,n and the independent coordinate x according to equations (ITU]) . (fTBD . 
flT^ . (ITS]) . (ITT]) . (ITT|) respectively) is determined by a solution s of ffTT]) . flT^ . We use 
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the following notation: Denote by Sa, Sjs etc the solntions of corresponding to 

parameters a, (3 etc. Denote the associated action on the helds by n —)■ Ua, m —)■ etc. 
Denote by SajS the solntion of flTTli . ffTSl) with u, m replaced by rria and parameter (3 (i.e. 
we start with a solution obtained from a BT with parameter a and are now considering 
acting upon it by a further BT with parameter f3). Denote the corresponding action on 
the helds by Ua Uai 3 , ma —>■ etc. 

The fundamental fact about double BTs, as proved in [60], is that they commute, i.e. 
Ua 0 = fnap = mpa etc. From, for example, the transformation law for the pCH held 
n, 031. it immediately follows that 


“1“ ^0.(3 ^(3 d* • 


(26) 


Checking the consitency of this with the versions of flTTD and flT^ satished by Sq,, s^, Sa/s, 
we obtain 


®a/3 ■Sa T 


{a - I3){al3 - SgSp) 
I3sq — asp 


SjSa — —SjS + 


(a - /3){a/3 - SgS/s) 
I3sa — (ysp 


(26) 


In fact it is possible to check directly that these formulas for Sa/s, s^a give solutions of the 
relevant versions of flTT]) and flT^ without any need to assume fl25l) . 

From fl26|) it follows that once and are known, it is possible to immediately hnd 
the action of a double BT. Using the transformation laws for m,p,n,x and fl26D we hnd 


Pa/3 


'^al3 


‘^0.(3 


^a/3 


af3 {{sg - Sfif - (a - /j)^) 
{(3sa - as^y 

{(3sa - aspY 
^ _ ^ (a - /3){a/3 - SgSp) 
I3sa — (ysp 


X + In 


Sp - Sa + a - /3 
Sp - Sa - a + 3 


For u and v we proceed as follows. From flTU]) and fllSp we obtain 


(27) 

(28) 

(29) 

(30) 


Ua + u - {Va — v)Sa = —2a (31) 

a 

and similarly 


U^+U + ^{Vl3-v)sp = 

1 

to 

(32) 

^al3 “h Uj3 d* 

—2a , 

(33) 

'^Q.j3 d“ Ua d" "pi^Ua^ Ua')Sa^ 

- 23 . 

(34) 


Eliminating Va^vp^Vap from these 4 relations, using fl26|) for Sap and ffTOj) for Ua,up we 
obtain 


2(a 

Ua/3=U - 


3) ((« 


3){<y + 3 3- u) + {sp — Sa){sp + Sa + Ux)) 
(a - 33 - {sp - Say 


(35) 
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(36) 


Similarly, by first eliminating the fields u, 


2(a 

Vap =V - 


/3) ((g - + {s^ - Sa)u + 2{asf} - ^Sg)) 

{a - I3y - [sfi - Sa)2 


Eqnations fl^ . fl^ . fl^ . fl30|) . fl35|) and fl36|l are algebraic formnlas for the implemen¬ 
tation of a donble BT given Sa and sp. However Sa and sp also determine the imple¬ 
mentation of the original single BTs, so it is natnral to try to eliminate them to obtain 
nonlinear snperposition formnlae for each of the qnantities p, m, n, x, u, v. For example, 
for X we have, from flT^ . 






a , 


Sfi 


e* -f 

gx _ g3;/3 


and nsing these in flHOD gives 

(e^ — a 

(e® — e^P){e^°‘ — 


(37) 


Thus we see e* satisfies the cross-ratio equation, equation Al[5 = 0] in the ABS classif- 
ciation [1]. Similarly for n we obtain 


fiifia + n — na){2a — -t- = aififi + n — np){2l3 — -|- nap) , (38) 


which is also the cross-ratio equation after a simple held redehntion. For p the situation 
is a little more complicated as we have 

4 = a(a + ppa) , sl = I3{I3 + ppp) , 

and knowledge of pa only determines Sa up to a sign. As a result, for given p,Pa,Pp 
there are 4 possibilities for pap, which are given by solutions of the two multiquadratic 
quad-graph equations 

4:a/3{a - /3){p - Pap){Pa - Pp) + afiip - Pap)‘^{Pa - Pp)‘^ 

+a(a- fi)(ppa-PpPap)"^ + fiifi - (y)iPPp - PaPap)"^ = 0, (39) 

-Aafiia + fi){p + Pap){Pa + Pp) - afiip + Pap)‘^{Pa + Pp)^ 

+a{a +/3){ppa - PpPap)'^ + fiifi + (y){ppp - PaPap)'^ = 0. (40) 


The hrst of these is precisely the H3* equation in the Atkinson-Nieszporski classihcation 
of integrable multiquadratic quad graph equations [2], as is the second after a simple held 
redehnition. 

For u and v we have not succeeded to write a single superposition principle not 
involving any of the other helds. However, using the relations flM]) - fl341) it is possible to 
write the following superposition principles involving, respectively, just u and n, and just 
V and n: 


a 


hua + u + 2 a Uap + up + 2 a\ ^ hup + u + 2/3 

\ Un, -n nr,ft - na ) ^ \ ur - n 




{vp - v){np - n) + {Vap - Va){naP - Ua) 




= 8(/5 — a) 


{Vg - v){na -n) + {Vgp - Vp){nap - np) 

a 

(42) 


Here the helds n satisfy the cross-ratio type equation fl38|) . 
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5 Travelling wave solutions 

In this section we apply the BT fllOI) . fll4p where s satishes dm),(Ell) to the constant 
solution of CH u = mq 7 ^ 0, to obtain travelling wave solutions, specihcally soliton and 
cuspon solutions. These and other travelling wave solutions have been extensively studied 
in the literature, see for example |10l HU EJ El Ell SZl SHI EEl EH [M], and the BT is just 
one of many methods to derive them. The advantages of the BT will become apparent 
when we study superposition in the next section. 

If a(Q! + Mo) > 0 there are two kinds of real solutions of flTT]) . f[T^ : 

r-r - ^( ^/a{uo + a){x-Xo + {2a - Mq) t) 

Sa = V ® i '^0 + «) tanh - 

\ ‘lOL 

which we call the “tanh-type” solution, and the same with tanh replaced by coth, which 
we call the “coth-type” solution. As we will see both of these give rise to travelling wave 
solutions. If a{a + Uq) < 0 then there are real solutions 

, - - -^—r ( A/-a (mo + a) (x - Xo + (2a - uq) t) 

Sa = V -OL (mo + «) tan --- 

\ lOL 




and the same with tan replaced by cot, and an overall minus sign. Both of these give rise 
to periodic solutions (see for example EEZ]), but these will not be studied here. 

Returning to the case a(a + Uq) > 0, it is useful to write a + Mq = af/^, where U > 0, 
so the solution flT5D becomes 

Sa = all tanh f ^ (x — Xq + (3 — f/^) at) j , (45) 


and the same with coth for a coth-type solution. Using flTOH . flTdD the resulting solution 
is Ua{xa, t) where 


Ur 


= a{U^-3) + 


2 a(f/2 - 1) 


Xrv = X 


In 


tanh^ (x — Xo + (3 — U^) at)) — 1 
U tanh (x — xo + (3 — U'^)at)) — 1 


U tanh (y (x — Xo + (3 — U‘^)at)) -|- 1 


(46) 

(47) 


or the same with coth. Finally, writing z = x — xq + [3 — at, the solution becomes 
Ua{xa,t) where 


u 


a 


Xq — Xo + (3 — U‘^)at 


a{U^ 


Z + 1x1 


I MU" - 1) 

^ U^iun\i^\Uz- 
U tanh ^Uz — 1 
U tanh z + 1 


1 ’ 


(48) 

(49) 


this being a tanh-type solution, or a coth-type solution, which is the same with tanh 
replaced by coth. Both tanh-type and coth-type solutions are travelling waves with speed 
c = (U^ — 3)a, written in an implicit form. The hrst step in analyzing these solutions is to 
decide whether the maps from Xa to z are bijections. For tanh-type solutions with U < 1, 
neither the factor in the numerator or in the denominator inside the In can vanish, and 

















thus Xa only tends to (plus or minus) infinity as 2 ; tends to (plus or minus) infinity. The 
corresponding solutions are solitons which tend to uq = a{U‘^ — 1) at spatial infinity, with 
speed c = a(17^ — 3) and central elevation —a(l + f/^) = c — 2uo. Note that since 


U = 




3^0 — c 
Mo - c 


and 0 < 17 < 1 we must either have c < 3 mo < 0 or 0 < 3 mo < c. Figure 1 displays 
the soliton profile for c = 2 and uq = 0.5, 0.1, 0.02. (For negative mq and c the soliton is 
inverted.) 



u0=0.5 - u0=0.1 u0=0.02 


Figure 1: Soliton profile, c = 2, wq = 0.5, 0.1, 0.02. 


Of particular interest is the limit of the soliton for fixed c and mq i 0 (for c > 0) or 
Wo t 0 (for c < 0). Figure 2 shows Xa as a function of 2 ; (for xq = t = 0), Wq, as a function 
of 2 ; and Ua as a function of x^ in the case c = 2, uq = 10“®. Xa is close to zero, and 
Ua is close to c for a large interval of 2 ; values of size O (|ln(no/c)|) around z = 0. In the 
plot of Ua against Xa this gives rise to a sharp peak. This is the peakon limit. To see this 
analytically it is possible to use fl^ to find z in terms of u (with a ± uncertainty as it 
is necessary to take a square root), and then fl4^ becomes 


Xa — xq + {3 — U‘^)at 



(50) 


Both terms on the RHS diverge as Uq —)■ 0, but it is straightforward to extract the 
divergent behavior, which cancels between the terms, and to obtain the limit, which is 
simply ±ln (^). 

Moving now to tanh-type solutions with 17 > 1, from fl4^ we expect Xa to diverge 
when iaxi\i]^Uz = and thus the map from z to Xa will not be a bijection. Figure 3 
shows Xa and Ua as functions of 2 ; for c = 2 and uq = 3. The map from 2 ; to Xa is 3 to 
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Figure 2: The soliton with c = 2, approaching the peakon limit. Uq = 10 x as a 
function of z, u as a function of 2 ; and n as a function of x. 




Figure 3: tanh-type solutions with c = 2 and uq = 3 {U = \/7 > 1), Xa and Ua as 
functions of z. The map from z to Xa is not 1 — 1. 

1 and thus there are 3 corresponding solutions of CH, depicted in Figure 4. Since these 
are all unbounded we do not devote further attention to them. 

Moving now to coth-type solutions, the situation is very similar, but now the map 
from 2 ; to Xa will be 1 — 1 if > 1 and many to 1 if f/ <1, and there is a subtlety 
arising due to the divergence of coth ^Uz at z = 0. For coth-type solutions with U < 1, 
Xa diverges when cot\i \Uz = ±^. The map from to Xa is once again 3 to 1. Figure 
5 shows Xa and Ua as functions of z for c = 2 and Uq = and Figure 6 shows the 
3 corresponding solutions of CH. The subtlety, as can be seen in Figure 7, is that the 
solution corresponding to the range of 2 ;’s that includes zero, has a cusp at z = 0, arising 
from the divergence of coth ^Uz. Since at this point Ua is not differentiable, it is necessary 
to ask in what sense this is a solution of CH. Fortunately, the value of Ua at the cusp is 
c, which makes it possible to interpret the solution in a weak sense [27], though we do 
not go into details here. 
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Figure 4: tanh-type solutions with c = 2 and Uo = 3 {U = \/7 > 1), as a function of 
Xa for the 3 unbounded solutions of CH. 




Figure 5: coth-type solutions with c = 2 and uq = ^ {U = ^ < 1), Xa and Ua as 
functions of z. The map from z to Xa is not 1 — 1. 

For coth-type solutions with 17 > 1, the map from z to Xa is a bijection, and once 
again there is a single solution of CH, but with a cusp at ^ = 0 — this is the cuspon 
solution. Due to the requirement 17 > 1 cuspon solutions only exist with speed c < uq 
if Uq is positive, and speed c> uo ii uq is negative. Figure 8 illustrates cuspon solutions 
with c = 2 for mq = “1; —0.5, —0.1. (For positive Uq the cuspon is inverted.) Note that 
the central elevation of the cusp is c, as required for it to be a weak solution. For c > 0 
(c < 0) it is possible to consider the limit of the cuspon as Mq t 0 {uq I 0), and this is 
once again the peakon limit. 
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Figure 6: coth-type solutions with c = 2 and Uq = {U = ^ < 1), as a function of 
Xa for the 3 unbounded solutions of CH. 



Figure 7: coth-type solutions with c = 2 and uq = \ {U = ^ < 1), close up on the cusp 
in one of the solutions of CH. 


We summarize the travelling waves presented in this section in the following table. 
All the solutions have asymptotic height uq-. 


tanh-type 

U < 1 

soliton 

central elevation c — 2uq 
Mo > 0: c > 3 mo 

Mo < 0: c < 3 mo inverted 


U > 1 

unphysical 


coth-type 

U < 1 

unphysical 



U < 1 

cuspon 

central elevation c 

Mo > 0: c < Uq inverted 

Mo < 0: C> Uq 
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uO=-l - u0=-0.5 - u0=-0.1 


Figure 8: Cuspon profile, c = 2, mq = —0.5, —0.1. 

6 Two wave solutions 

The first investigations of two wave solutions were [SU] and |22], both of which required 
some element of numerical computation. However, since then, a substantial literature 
[3ll|39lll9l[l8l[Mlll3l[Sllia[52l[l9l[63] has developed on multisoliton, multicuspon and 
soliton-cuspon solutions. The known methods for analytic construction of solutions in¬ 
clude a determinantal formula based on the inverse scattering approach, a Hirota bilinear 
form for CH and a reciprocal transformation relating the CH hierarchy to the KdV hierar¬ 
chy. (For multipeakon solutions very different techniques are involved [9l HI El [5], [51].) As 
we will shortly see, use of the superposition principle gives a further very simple method. 

In our approach, two wave solutions should be obtained using formulas fi^ and flHUD . 
taking u = Uq to be constant and Sq (s^) either of tanh-type, as given in (H5]) where 

U = Ua = {U = Ufj = and xq = (a^o = or of coth-type, 

which is identical but with coth. The only question is which superpositions of this type 
give maps from x to Xqis that are 1-1. 

Proposition. The following 3 superpositions give maps from x to x^p which are 1-1: 

1. tanh-type solutions Sa with [/„ < 1 with tanh-type solutions S/3 with > 1 (so 
^ < 0 < ^) — soliton-cuspon superpositions. 

2. tanh-type solutions with Ua < t with coth-type solutions with Up < 1, with 
Ua < Up — soliton-soliton superpositions. 

3. tanh-type solutions Sa with Hq > 1 with coth-type solutions sp with Up > 1, with 
Up < Ua — cuspon-cuspon superpositions. 

Note here, for example, that a soliton-soliton superposition is not as we might expect, the 
superposition of two tanh-type solutions with 17 < 1, but the superposition of a tanh-type 
solution with 17 < 1 with a unphysical coth-type solution with 17 < 1. 
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Proof. It is necessary to show in each case that neither the numerator or denominator 
of the expression inside the In in flHOj) vanishes, i.e. that |s^ — So,| 7 ^ |/3 — Q!|. In the 
calculations below we repeatedly use the identities 


a = 


Uq 


ui-i^ 


/9 = 


Mo 


Uj-l 


1 . In this case we have 


|s /3 — Sq,| = tanh (...) —of/a tanh (...) 

< + \(y.Ua\ as I tanh I < 1 


mo 17/3 

1 

UoUa 

C/|-l 




= kol 
= l^ol 
< kol 


Up 




U }-1 Ul -1 

1 1 

+ 


as 0 < I/q, < 1 < 17/3 
1 1 


17 | - 1 f //3 + 1 f /2 - 1 f/„ + 1 

1 1 
as —-- > 


Uj-I Ul-I 


Ua + 1 Ufj + 1 


= \(3 — a\ . 


2 . In this case we have 

kal + I/S — a| = |al 7 a tanh (...) I + 1/3 — a| 

< \aUa\ + 1/3 — aj as I tanh | < 1 


UoUa 

+ 

Mo 

Mo 

ui-i 

Uj-l 

ui-i 


i«„i I ^ 


1 - C/J 1 - £/| l-U^ 


as 0 < f/„ < 17/3 < 1 


= |wo| 


1-172 1 + 17 , 


|mo| I . + 


1-172 ' l + f //3 1 + f/, 


< |mo| 


Up 


1 1 
as --— < 


1-172 1 + 17/3 1 + 17 , 


= 

< 1 / 317/3 coth (...) I as I coth | > 1 

= k/3l • 


This contradicts |s /3 — Sa| = |/3 — a|, as the latter implies |s/ 3 | < |sq| + |/3 — a|. 
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3. Similarly to case 2 we have 


Sal + \(3 - a\ = 
< 


< 


< 


\aUa tanh (...) | + |/9 — a| 

\aUa\ + I/S — a| as I tanh | < 1 


UqU a 

+ 

Uq 

Uq 

f /2 _ 1 

^ OL 

S' 

1 

172-1 

^ OL 


+ (72-l) 


5:tt) 

^ J Ub 1 1 A 

V^|-1 Uf, + i^ u^ + ij 

11^/3 1 1 

\(3Uj3 coth (...)! cLs| coth I > 1 

k/?! • 


as 1 < 17/3 < 17 q 


It remains to present the plots of some superpositions. Figure 9 shows a tanh-tanh 
superposition with uq = 1, Xa,o = 0, X/3,0 = 2, Ca = uq — 2a = 4, cp = uq — 2f3 = —1. 
(For mq > 0 such solutions exist provded Cq > Smq and cp < Uq — note here that cp 
can be positive or negative.) Figure 10 shows a tanh-coth soliton-soliton superposition 
with Uq = 1, Xa,o = 0, xpfi = —10, Ca = Uq — 2a = 4, Cp = Uq — 2(3 = 6. Figure 
11 shows a tanh-coth cuspon-cuspon superposition with uq = 1, Xa,Q = 0, xp^ = —2, 
Ca = Uq — 2a = —1, Cp = Uq — 2(3 = —2. 


^=-4, -1,4000 ^ = 0.40000 t = 4,6000 


A ^ 
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J\ 



( 

1 
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-1 
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Figure 9: Tanh-Tanh Superposition, uq = 1, Xa,Q = 0, xp^ = 2, Ca = 4, cp = —1. Plots 
for times t = —4, —1.4, 0.4,4.6 from left to right. 


Note that in all the plots we have taken uq > 0, in which case, as in the previous 
section, the soliton solutions have positive speed c > 3 mo and central elevation c — 2uq > 
Uq, whereas the cuspon solutions have speed c < uq, which can be positive or negative, 
and central elevation c < Uq (i.e. they might be called “anticuspons”). In the case Mq < 0 
everything is inverted and reversed (reflecting the t —)■ —t, u —?■ —u symmetry of ([2]): 
solitons have negative speed c < 3 mo and central elevation c — 2uq < Uq (antisolitons), 
and cuspons have speed c> uq and central elevation c> uq. 
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Figure 10: Tanh-Coth Soliton-Soliton Superposition, uq = 1, Xa^ = 0, = —10, 

Ca = 4, C/3 = 6. Plots for times t = —6, —1, 5, 8 from left to right. 





Figure 11: Tanh-Coth Cuspon-Cuspon Superposition, uq = 1, Xa,o = 0, = ~2, 

Ca = —1, = —2. Plots for times t = —10, —5, —2,4 from left to right. 

7 Symmetries and conservation laws for the Camassa- 
Holm equation 

In this section we show how to use the BT to obtain inhnite hierarchies of symmetries 
and conservation laws for CH and pCH, following the general methodology described in 
|55j . The discussion of symmetries of CH in the literature is limited, though the existence 
of an inhnite number of symmetries is implicit from the bihamiltonian structure given in 
[9] . In the series of papers [56l EU EH EH [27] , Reyes and collaborators present nonlocal 
symmetries of CH depending on a parameter, and then expand in powers of the parameter 
to obtain local symmetries, though limited details are given. Some explicit formulae 
appear in [21]. Our approach is related, but we will not discuss the connection explicitly. 

As a starting point for our discussion of symmetries we could take the generating 
symmetry ([9]) for aCH, and work out the induced action on x, the independent variable 
in CH. But a more direct approach is to look at the superposition principle fl35|) . fl30|) in 
the limit that {3 tends to a, but Sjs tends to a second solution of flTTD - flT^ distinct from 
Sa- More explicity, setting = a — |e, Sq, = Sa\ = Sa^ -1- 0(e) in we obtain 




M -1- e- 


X -f- 


Sq ^ + Sa^ + Ua 

d2) _ „(1) 

OQ; Oq 


+ 0{e^) , 


A2) _ Ab 

O/y Ory 


0(e^ 
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We deduce the generating symmetry for CH X = where 


Q" = 


.(2) _ Jl) ’ 

oq '-’a 


Q“ = 


sL ^ + sL ^ + Ma; 
J2) _ Jl) 

Oq: '^Q! 


(51) 


Here Sa\ Sa^ are two different solutions of ffTTli . ffT^ for the same parameter a. This 
symmetry depends upon a; expansion in a (formal) power series in a will give an inhnite 
hierarchy of symmetries. However before we do this, we exploit the fact that a gener¬ 
alized symmetry of the form X = which acts on both the dependent and 

independent variables can be transformed to a generalized symmetry which acts only on 
the dependent variable |16] with characteristic Q = — Q^u^- Here we have 


Q 


Sa ^ + So * 




(52) 


This is also the characteristic for a symmetry of the full family of equations flT^ . 

The next thing to do is to hnd a (formal) asymptotic series solution of flTT]) - flT^ for 
small |a|. This takes the form 

OO 

Sa = '^Sna^, (53) 

n=l 

where 


Si = y/m, S2 


^1,X 


^n,x 

Si 


n—2 

^ j ^n,2 ^ ^ Si-\-2Sn—i 

i=0 


n = 2, 3.... 


A second solution of ffTT|) - ffT^ can be obtained by replacing by —a^. So we get 


OO 

n=l 


n=l 


Plugging this into we obtain 


Q _ E^=iS2na” 

/— 

1 S2n—1^ 


(54) 


(56) 


The expansion of fl5^ around a = 0 gives an inhnite hierarchy of symmetries of CH. The 
hrst few of these take the form 


X2 

^3 



^mmxx 


— 5ml + 4m2 
-^2 


d 


du 


(56) 

(57) 


/ QAm^rrixxxx — 4:4:8m‘^mxmxxx — IQOm^mxx + \ 

-336m2mL + 280m2m2 - 1155m^ - 48m^ 


du 


(58) 


V 
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The fact that all the characteristics are x-derivatives is indicative that these symmetries 
can be derived from corresponding symmetries of pCH. The generating symmetry for 
pCH (up to an irrelevant overall constant factor) is thus where 




1 




(59) 


As stressed before, the symmetry with characteristic 0521) is a symmetry for the full 
family of equations (IT^ . including the HS equation. For HS the asymptotic series solu¬ 
tions of fl23l) - fl^ takes the form 


OO 

Sa ^ ^ 

n=l 


where 


• ! - 

5l ^2 : '^n+1 


= (60) 

n=l 



Proceeding as before gives an inhnite hierarchy of symmetries for HS, with the hrst few 
taking the form 


^2 

^3 



^XX 


^XX 


-f448 


'^xxx'^xxxxx 

9/2 

Uxx 


+ 336 


ut 


9/2 

Uxx 


64 


Ur 


7/2 

'Uxx 


X 


A 

du 


Using the fact that if a single solution of the Riccati equation flTT]) is known then it is 
possible to hnd the general solution by quadratures, it is possible to rewrite in the 
form 


Q — 1 + 

and fl55]) in the form 




a 




• XQ 


Q'^ = I + (Ue“ . 


1 


fXQ 


Here C is an arbitrary constant. Since a linear combination of symmetries is a symmetry, 
both terms on the RHS are by themselves the characteristics of symmetries. The sym¬ 
metry associated with the factor multiplying C is the nonlocal symmetry hrst presented 
in [56]. The relation between Backlund transformations and nonlocal symmetries has 
recently been discussed in 021 . 

A conservation law (CL) for a PDF for the scalar function u{x,t) is an expression 


Tt + Xx — 0 


which holds on solutions of the equation. Conservation laws for CH can be obtained from 
flT^ by writing it in the form 

St + {su - a{ux + 2s))^ = 0 (61) 
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Using the expansion fl5^ for s in fl6T]l we obtain an infinite hierarchy of conservation laws. 
Terms with integer powers of a in this expansion give trivial CLs. To prove this, observe 
from (15T|) that the terms with integer powers are obtained by setting s = | ^ + Sa ^ j 

in (I6T]) . Bnt from flTT]) it is simple to verify that 

^ (4^^ + 4^^) = -«(In 14^^ - 4^4x • 

Thus to obtain nontrivial laws we look at only the half-integer powers of a. Thus we set 
s = — Sq ^ in fIM]) to obtain the generating conservation law 

T = 4"-4"'. (62) 

X = (ti-2a)(4‘’- 4’’) ■ (63) 

The expansion of the generating conservation law around a = 0 gives an inhnite hierarchy 
of nontrivial CLs for CH and pCH. The hrst few take the form 


Ti = \/m , 

Xi = uTi , 

^2 = + 4mm,, - 5m^) , 

X2 = UT2 — 2 Ti , 

T 3 = ]^28m^V2 (64m^m,,,, - llOSm^ 1768mm^m,, - 304m^m^^ - 448m^m,m,,, 

—96m^m,, -|- 200m^m^ — 16m^) , 

X3 = MT3 — 2T2 . 

Similar results can be obtained for HS and the full family of equations (ITU]) . The existence 
of an inhnite number of conservation laws for CH follows from the bihamiltonian structure 
for CH discovered in [9]. The local form of these conservations laws was hrst obtained 
in [23], and they were subsequently further studied, and their derivation simplihed, in 
numerous works such as [SHl[SZl[2Sl[2SlEZlE51[IIl[3ni[21|. The derivation given here can 
be easily related to previous ones, though the use of 2 solutions of flTTD - ffT^ to understand 
the triviality of “half” of the conservation laws, is, we believe, new. 

8 Concluding remarks 

In this paper we have explored the theory of the Backlund transformation for the Camassa- 
Holm equation. This is an unfamiliar type of BT, as it acts on one of the independent 
variables, as well as the dependent variables. However, it has emerged that it is just as 
useful — using the superposition principles for the action on the different variables, we 
can exploit the BT to write down two wave solutions, just as is done for standard inte- 
grable equations such as KdV. Furthermore, we have shown how a double BT encodes an 
inhnite set of symmetries for CH, and the relationship of the BT and conservation laws. 

We have seen that the BT can also generate “unphysical” solutions, by which we 
mean solutions for which the new independent variable is not a 1 — 1 function of the 
old independent variable. Going beyond two wave solutions, it is not clear exactly what 
superpositions are allowed without creating singularities, though it seems to be a rea¬ 
sonable hypothesis that all possible combinations of solitons and cuspons can be formed. 
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with the speeds permitted by the value of mq, as listed in the table at the end of section 
5. It seems to us that this is a problem that remains to be handled indepdendent of the 
method used for constructing multiwave solutions. 

Peakons emerge from both solitons and cuspons in the limit Mq 0 (with one giving 
rise to peakons of positive speed and one to peakons of negative speed, depending on 
whether the limit is taken from below or from above). This is an extremely singular limit. 
We have not yet found a way to apply a superposition principle directly to peakons, but 
we continue to search. 

Finally, one more general comment. The BT, in its minimalist form, is the transfor¬ 
mations (ini]) and (ITT]) where s satishes flTT]) - ffT^ . The latter equations for s are equivalent 
to the Lax pair, or linear system, for CH. So the BT seems to be more than the linear 
system. We wonder if there is a case of an integrable system without a BT? 
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